Sensitivity impacts owing to the variations in the type of zero-range pairing forces on the fission 
properties using the density functional theory* 
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Using the Skyrme density functional theory (DFT), potential energy surfaces (PES) of 7“°Pu with constraints 
on the axial quadrupole and octupole deformations(g20 and q30) were calculated. The volume-like and surface- 
like pairing forces, as well as a combination of these two forces were used for the Hartree-Fock-Bogoliubov 
(HFB) approximation. Variations in the least-energy fission path, fission barrier, pairing energy, total kinetic 
energy, scission line, and mass distribution of the fission fragments based on the different forms of the pairing 
forces were analyzed and discussed. The fission dynamics were studied based on the time-dependent generator 
coordinate method (TDGCM) plus the Gaussian overlap approximation (GOA). The results demonstrated a 
sensitivity of the mass and charge distributions of the fission fragments on the form of the pairing force. Based 
on the investigation of the neutron-induced fission of 2°°Pu, among the volume, mixed, and surface pairing 
forces, the mixed pairing force presented a good reproduction of the experimental data. 
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I. INTRODUCTION 


Nuclear fission, i.e., the phenomenon that one (usually 
heavy) atomic nucleus may separate into two or more frag- 
ments, has been discovered for more than eighty years [1, 2]. 
It is accompanied by the release of abundant energy [3] and 
has a wide range of applications. In addition to the important 
applications in energy and production of rare isotopes, fission 
also plays a crucial role in fundamental physics, such as syn- 
thesizing superheavy elements [4—6], as well as constraints on 
the r-process in neutron-star mergers [7] et al. Therefore, sev- 
eral theoretical approaches have been utilized to describe the 
fission process and observations [8—17]. As shown in Ref. [8], 
the microscopic models that were applied to fission thus far 
utilize the density functional theory (DFT), which is based on 
effective nucleon — nucleon interactions. 

As one of the dominant residual correlations in the atomic 
nuclei, the pairing interaction is critical for understanding the 
fission process. Extensive studies regarding the influence of 
pairing interactions on the fission properties were performed, 
such as the effect on the fission barrier heights, fission iso- 
mer excitation energies, and collective inertia [18-23]. In 
Ref. [20], the fission dynamic calculation based on the co- 
variant DFT was performed for the fission of 2?°Th, in which 
both the symmetric and asymmetric fission modes co-exist. 
The asymmetric fission mode dominated as the pairing force 
decreased, whereas the symmetric fission mode dominated as 
the pairing force increased. The time-dependent superfluid 
local density approximation (TDSLDA) method has demon- 
strated that fission significantly accelerates as the pairing 
force increases [24]. The effect of the dynamic pairing corre- 
lations on the fission process was studied in Refs. [25, 26]. 
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In our previous work, the role of the pairing force on the 
static fission properties and fission dynamic process were 
studied [22, 27-29], which demonstrated that the variation in 
the strength of the pairing correlation can significantly influ- 
ence the fission process. However, whether the form of pair- 
ing force can also impact the fission properties is noteworthy. 
A schematic pairing force was originally introduced in [30], 
and parameters of the density-dependent pairing correlation 
were studied in Refs. [31-33]. Specifically, considering 
Skyrme-DFT, different types of pairing forces have been used 
for studying the nuclear structure, such as the volume-like, 
surface-like, and mixed-type pairing forces. In Refs. [34-36], 
various types of the pairing interactions (volume-, surface-, 
and mixed-type pairing forces) were used to study the pair- 
ing gaps in even-even nuclei over the entire nuclear chart, the 
odd-even staggering behavior of binding energies around tin 
isotopes, and to predict the two-neutron separation energies 
and neutron pairing gaps[36]. In this study, whether the form 
of the pairing force can influence the static aspects of the fis- 
sion properties and dynamics was investigated. The fission 
process of 74°Pu was studied in Skyme DFT and the time- 
dependent generator coordinate method (TDGCM) + Gaus- 
sian overlap approximation (GOA) (TDGCM + GOA) frame- 
work. 

In Section II, we briefly describe the main features of the 
theoretical approach. Details regarding the calculated results 
for the least-energy fission path, fission barrier, pairing en- 
ergy, total kinetic energy, scission line, and mass distribu- 
tion of fission fragments using various types of pairing forces 
are analyzed and discussed in Section II. Finally, Section IV 
presents a summary of the principal results. 


Il. THEORETICAL FRAMEWORK 


Skyrme DFT was applied as the microscopic method to 
study the static fission properties and prepare the input files 
for the dynamic calculation. The dynamic process was fur- 
ther investigated in the framework of TDGCM. Thus, in this 


section, these two methods are briefly explained. Detailed 
formulations of the Skyrme DFT can be found in Ref. [37], 
and those of TDGCM can be found in Refs. [38-40]. 


A. Density functional theory approach for the description of 
PES 


For the local density approximation of DFT, the total en- 
ergy of the finite nuclei can be calculated using the space in- 
tegral of the Hamiltonian density 1(r), which consists of the 
kinetic energy 7, potential energy x+, and pairing energy X+ 
densities: 


H(r) = x(n) + ut UO). O 


Here, 7(r) is the density of the kinetic energy, and the 
symbol t = 0,1 indicates the isoscalar or isovector, respec- 
tively [41]. 

The mean-field potential energy in the Skyrme DFT usu- 
ally has the following form: 
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Here, the particle density p+, kinetic density 7, and spin cur- 
rent vector densities J;(t = 0,1) can be obtained from the 
density matrix p;(ro,r’o’), depending on the spatial (r) and 
spin (g) coordinates. In the aforementioned formula, C? T 
C£”, and etc. are the coupling constants in the Hamiltonian 
density H(r), corresponding to the different types of densi- 
ties, most of which are real numbers; C?? = Ch? + CPF od 
is an exception, which is the traditional density-dependence 
term. Expressions relating the coupling constants to the stan- 
dard Skyrme parameters can be found in Ref. [42]. Specif- 
ically, the spin-orbit interaction in the Skyrme force corre- 
sponds to the term CPY 7 p,V - Jy. 

In the DFT, the pairing correlation is usually incorporated 
by the Hartree-Fock-Bogoliubov (HFB) method [37]. For the 
Skyrme energy density functional, a commonly used pairing 
force is the density-dependent, zero range potential, which 
can be expressed as follows: Refs. [8, 43]: 


Vair, 1) z Vier) [1 ma yY] -r), (3) 


Here, ven) is the pairing strength for neutrons (n) and pro- 
tons (p), the exponent y of the density-dependence affects the 
appearance of the neutron skins and halos [44], for which 
y =1 is widely used [45, 46]. po is the average density in- 
side the nucleus (often considered as the saturation density of 
nuclear matter, and set as 0.16 fm~3), and p(r) is the total 
density. Different types of pairing forces can be obtained by 
choosing different values of 7. The pairing force is ” volume- 
like”, which indicates that if 7 = 0, there is no explicit den- 
sity dependence. The pairing force acts equivalently inside 
the nuclear volume. In contrast, the pairing force will be 


“surface-like” for 7 = 1, which has a significant effect around 
the nuclear surface and a small impact in the center area of 
the nucleus. The choice of 7 = 4 is often called a mixed- 
pairing force, which is the average of these two types of pair- 
ing forces. To test the sensitivities on the fission-related prop- 
erties based on the form of the pairing force, we considered 7 
= 0, 0.25, 0.5, 0.75, and 1. 7 = 0.25 and 0.75 were used for 
test purposes, and the other choices of the 7 values have been 
frequently used for structural studies. 


B. Time-dependent generator coordinate method for fission 
dynamics 


Fission is a large-amplitude collective motion and can be 
described as a slow adiabatic process driven by a few collec- 
tive degrees of freedom within the framework of TDGCM. 
In this approach, the many-body state wave function of the 
fissioning system can be expressed by the following generic 
form: 


w(t) = I f(a,t)|®(q))qa. (4) 


where |®(q)) consists of a set of known many-body wave 
functions parameterized by a vector of continuous variables 
q. Each of these q is a collective variable and is cho- 
sen based on the specific physics problem. The quadrupole 
moment Cais and octupole moment Oana are usually chosen 
as the collective variables for the fission study. f(q,t) is 
the weighted function, which is solved by using the time- 
dependent Schrodinger-like equation in the space of the co- 
ordinates q. Under the GOA, this equation can be expressed 
as follows: 
og (q, t) a 


iñ = Aon (q) (4, t). (5) 


The collective Hamiltonian Heo (q) is as follows: 
f h? ð ð 
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where V (q) is the collective potential, and the inertia tensor 
Bi;(q) = M~1(q) is the inverse of the mass tensor M. The 
potential and mass tensor are determined by the Skyrme DFT 
in this study. g(q, t) is the complex collective wave function 
of the collective variables q and contains all the information 
regarding the dynamic of the fission system. 

For the description of fission, the collective space is divided 
into an inner region with a single nuclear density distribution, 
and an external region that contains the two fission fragments. 
The scission contour defines the hyper-surface that separates 
the two regions. The flux of the probability current through 
this hyper-surface provides a measure of the probability of 
observing a given pair of fragments at time t. The integrated 
flux F'(€,t) for the surface element € on the scission hyper- 
surface is calculated as follows: 


t 
F(é,t) = f it f Nadas, (7) 
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as in Ref. [40], where J(q, t) is the current 


h * * 
Ja t) = 5-Bla)[s" (a. 4) Vela, t) -ela t)Ve" (a, t)]. (8) 
The yield of the fission fragment with mass number A can be 
calculated as follows: 


Y(A)=C) lim F€,t) (9) 
EEA 


where A indicates a set of all the surface elements € on the 
scission hyper-surface with a fragment mass of A, where C 
is the normalization constant to confirm that the total yield 
is normalized to 200 as usual. The yield of the fission frag- 
ment with the charge number Z can also be obtained from 
the integrated flux F'(€,t). The mass number A can be re- 
placed with the charge number Z in Eq. 9; the summation is 
then over the set of all the € values on the scission frontier 
with the fragment charge number Z. In this study, the FE- 
LIX(version 2.0) [47] computer codewas used for modeling 
the time evolution of the fissioning nucleus in the framework 
of TDGCM+GOA. 


I. RESULTS AND DISCUSSION 


The influence of the different types of pairing forces 
on the fission properties were studied based on Skyrme 
DFT with SkM* parameters [48]. The DFT solvers HF- 
BTHO(V3.00) [49] were used to calculate the potential en- 
ergy surfaces (PESs). An axial symmetry was assumed. 
Thirty-one major shells of the axial harmonic-oscillator 
single-particle basis were used, and the number of the ba- 
sis states were further truncated to be 1100. In this study, we 
considered five values of 7 (0, 0.25, 0.5, 0.75, and 1) in Eq. 3 
as the different types of pairing forces. For each value of 7, 
the pairing strength for the neutrons and protons was adjusted 
to reproduce a pairing gap of 74°Pu extracted from the three- 
point formula of the odd-even mass staggering. A cutoff of 60 
MeV was used as the pairing window in all the calculations. 

Fig. 1 presents the pairing strengths of the neutrons and 
protons for the pairing forces with different values of 7. For 
7 near 0, the pairing tends to occur equivalently in the nu- 
clear volume. When it is near 1, the pairing tends to peak at 
the nuclear surface. For 7 between the values of 0.0 and 0.5, 
the absolute value of the pairing strength increased nearly lin- 
early. However, for 7 = 1.00, there was a sudden increase in 
the pairing strength. The surface-like pairing force requires a 
significantly larger strength to produce the same pairing gap 
compared to the volume-like or mixed-type pairing forces. 
These results are consistent with previous studies. As demon- 
strated in Refs. [50, 51], the surface pairing force is also sig- 
nificantly stronger than the volume and mixed pairing forces. 


A. Potential energy surface 


As indicated in Sec.ILB, considering the adiabatic approx- 
imation approach for fission dynamics, obtaining the precise 
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Fig. 1. (color online) Pairing strength of the neutron(square) and 
proton (circle) as a function of the parameter 77. The strength Vo was 

adjusted to reproduce the empirical pairing gaps in 74°Pu. 


multidimensional PESs is the first step toward the dynamical 
description of fission. In this study, we chose the quadruple 
moment (q20) and octuple moment (q30) as the collective pa- 
rameters, which are the most important collective degrees of 
freedom for the nuclear fission study; they describe the elon- 
gation of the nucleus and mass asymmetry, respectively. Fig- 
ure 2 presents the PESs of ?4°Pu calculated by using the HFB 
method with five different types of pairing forces (7 = 0, 0.25, 
0.5, 0.75, and 1.0) in the collective space of (q20, q30). The 
collective variables ranged from 0 to 600 b for goo, and from 
0 r 60 b?/? for qso with the step of Aq20 = 2 b and Aq3o=2 
be/?, 

As shown in Fig. 2, there is no notable difference in the 
topological properties of PESs with different types of pair- 
ing forces. Double-humped fission barriers were predicted 
for all the cases. An inner symmetric fission barrier followed 
by an outer asymmetric barrier was clearly distinguished. At 
q20 >200 b, symmetric valleys with large elongations were 
found. The symmetric and asymmetric fission valleys were 
well-separated by a ridge from (q20, 930) ~ (150 b, 0 b?/?) 
to (350 b, 20 b®/?), and the height of the ridge gradually 
decreased as the 7) value increased. Therefore, the density- 
dependent surface pairing force led to the reduction of the 
ridge height. In addition, the asymmetric fission channel was 
favored for all the least-energy fission pathways, as indicated 
by the red lines in Fig. 2. 

Energies of the symmetric and asymmetric fission paths 
as a function of the quadruple moment (q20) are provided in 
Figs. 3(a) and (b), respectively. The value of 7 varied from 
0 to 1, indicating that the *volume-like” pairing force transi- 
tioned into a surface pairing force. Fig. 3 demonstrates that 
the fission barrier heights and isomeric-state energy decrease 
as ņ increases. Specifically, when 7 > 0.5, the fission barri- 
ers explicitly decrease. For the least-energy fission pathway 
shown in Fig. 3(b), a smaller quadruple moment is needed for 
the occurrence of the scission for a larger 7. 

Table 1 lists the energies of the ground state, isomeric 
states, and fission barrier heights for the different types of 
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Fig. 2. (color online) The PESs of 240Pu in the (q20, q20) plane 
calculated with the SkM*-DFT and the five different types of pairing 
forces in the HFB approximation. The value of 7 for the different 
types of pairing forces is provided in each panel. The energies are 
provided relative to the ground-state energy. The contours join the 
points on the surface with the same energy, and the interval of the 
contour is 1.0 MeV. The least-energy fission pathway is indicated by 
a solid red curve. 


pairing forces, along with the corresponding quadruple and 
octuple moments for each state. The energies of the isomeric 
state and heights of the fission barrier decrease as the 7) value 
increases. Owing to the lack of triaxial deformation in the 
collective space in our calculation, the heights of the fission 
barriers would be higher than those experimentally obtained, 
especially for the inner fission barrier [52, 53]. As indicated 
in the table, for 7=1, the inner fission barrier height was near 
that demonstrated by the data, and the outer fission barrier 
was lower than that in the data, leaving no room for the triax- 
ial degree of freedom. Thus, 7=1, that is, the surface pairing 
force, may not be a good choice for the fission study. Based 
on Table 1, the deformations of these states, including the 
ground state, isomeric state, and inner and outer barriers, are 
generally not influenced by the type of pairing force. These 
deformations are mainly determined by the shell structure 
given by the mean-field potential. In our previous study [22], 
we also found that these deformations were relatively stable 
against the variations of the pairing strength. 
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Fig. 3. (color online) Energies along the symmetric and least-energy 
fission pathways in °*°Pu, with different types of pairing correla- 
tions, are shown in panels (a) and (b), respectively. All these ener- 
gies are relative to their ground-state values. 


B. Pairing energies 


Fig. 4 presents the pairing energies at different deforma- 
tions for various types of pairing forces. The pairing energies 
at the ground and isomeric states are smaller than those at 
the fission barriers. At the same state, the pairing energy in- 
creases as the value of 7) increases, especially for n=1, which 
is the surface-type pairing force. The pairing gaps at different 
deformations are provided in Fig. 5. Once again, the pairing 
gap has a minimum at the ground state and a second minimum 
at the isomeric states. The pairing gaps are large around the 
fission barriers. For the adjustment of the strength of the dif- 
ferent pairing forces used in this study, the same values of the 
pairing gaps in ?4°Pu were used. The figure demonstrates that 
at the smaller deformation region, the pairing gaps from the 
different types of pairing forces were relatively similar. How- 
ever, when the deformation was large, explicit discrepancies 
appeared (q20 >150 b). For the pairing force with a smaller 
7 value, the neutron and proton pairing gaps were generally 
smaller. 


C. Mass tensor 


The mass tensor M reflects the response of the fissioning 
system to the collective coordinate changes. In this study, 
mass tensors were obtained from the static calculations by 
using the Skyrme DFT with the perturbative cranking approx- 
imation. As shown in Fig. 6, the elements of the mass tensor 
Mo2, M33, and M3 were plotted as functions of q2ọ along the 
lowest-energy fission path for the different choices of 7. The 


TABLE 1. Energies of the ground state, isomeric state, and fission barrier heights. The empirical values are indicated by “*”. The quadruple 
and octuple moments at the corresponding states are provided. The labels 0.0, 0.25, 0.5, 0.75, and 1.0 indicate different types of pairing forces 


and are the values of 77. 


Ground state 


Isomeric state 


Inner barrier Outer barrier 


Exp -1813.45 [54] 2.8 [55] 6.05* [56] 5.15* [56] 
0.0 -1805.06 2.83 9.83 6.90 
Energy/MeV 0.25 -1805.15 2.63 9.58 6.64 
0.5 -1805.35 2.59 9.32 6.46 
0.75 -1805.68 2.31 8.72 5.84 
1.0 -1806.44 1.55 7.12 4.36 
0.0 (30,0) (86,0) (54,0) (124,8) 
Deformation 0.25 (30,0) (86,0) (54,0) (124,8) 
(q20/b,q30/b3/2) 0.5 (30,0) (86,0) (54,0) (126,8) 
0.75 (30,0) (86,0) (54,0) (124,8) 
1.0 (80,0) (86,0) (54,0) (124,8) 
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Fig. 4. (color online) Pairing energy of the ground state, fission iso- 0.00 


mer state, and inner and outer fission barriers in 240Pu as a function 
of the 7 parameter in the pairing force. 
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Fig. 5. (color online) Pairing gap along the least-energy fission path 


for the neutron (a) and proton (b) with different types of pairing 
forces. 


O 50 100 150 200 250 300 O 50 100 150 200 250 300 350 
,9(b) q,,(b) 

Fig. 6. (Color online)M22 (top panel), M23 (middle panel), and 

M33 (bottom panel) of the mass tensors are shown as the elongated 

deformations along the static fission path for different 7 values. The 

resulst in the left column are obtained by using GCM method, and 

those in the right column are derived by the ATDHFB method. 


mass tensor given by ATDHFB was larger than that by GCM 
for all types of pairing forces. As indicated in Ref. [8], this 
was caused by the missing correlations in the GCM method. 
Generally, as the quadruple moments increase, M22 and M33 
gradually decrease. Mo3 is negative, and its absolute value 
increases and has explicit fluctuations when the deformation 
is large. For different types of pairing forces tuned by 1, the 
mass tensor apparently decreases and demonstrates reduced 
fluctuation against deformations when 77 is large. As indi- 
cated in Ref. [28], considering the mixed pairing force of 
n = 0.5, the mass tensor decreases and fluctuates less when 
the pairing strength is decreased. The systematic behaviour 
of the 7 > 0.5 pairing force resembles that of increasing the 
strength of the pairing force. 
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Fig. 7. (color online) The scission lines using the criterion of gv = 
4 are shown for the different 7) values. 
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D. Scission lines 


Determining the scission frontier is critical for describing 
the fission dynamics. In DFT, the operator gy = (Qn) = 


(ew) is often used to evaluate the neck size of the fis- 
sioning nuclei. ay = 1 fm is often chosen. zy is the neck 
position, which has the lowest density between the two frag- 
ments. Generally, the neck size smoothly decreases as the fis- 
sioning nucleus elongates, and decreases to nearly zero after 
the scission, where the two fragments are sufficiently sepa- 
rated. In this study, we chose qy = 4 as the critical value for 
determining the scission line in 74°Pu, which has been used 
for 74°Pu in Refs. [28, 40]. This value was chosen at the edge 
of the sudden decrease in the neck size, maintaining most 
of the pre-fission configurations for further fission dynamic 
calculations. The scission lines in the PES of the (q20, q30) 
collective spaces obtained by DFT using the different types 
of pairing forces are shown in Fig. 7. Generally, the scission 
contours of the different 7 values display similar patterns. For 
different 7 values, in or around the region of the symmetric 
fission, the increase of 7 leads to a smaller quadrupole mo- 
ment at the scission point. Considering 7 = 1.0, symmetric 
fission occurred at the quadrupole moment q20 ~ 480 b for 7 
=1.0, and at approximately 550 b for the other 7 values. The 
shortest elongation occurred at q20 ~ 300 b when the asym- 
metry increased to the octupole moment q39 ~ 30 b?/?. Sub- 
sequently, the scission lines turned toward the upper-right di- 
rection until significant asymmetry appeared. The pre-fission 
region for 7 = 1 is explicitly smaller than the other cases. 


E. Total kinetic energy 


An important quantity in induced fission is the total ki- 
netic energy (TKE), which is obtained by the fission frag- 
ments. In this study, the total kinetic energy of the two sep- 
arated fragments at the scission point can be estimated as the 
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Fig. 8. (color online) Total kinetic energies of the nascent fission 
fragments as functions of the heavy fragment mass for 74°Pu based 
on the different 7 values. Open symbols represent the calculated 
results, and solid symbols indicate the experimental data [61—63] 
for the thermal neutron. 


dicates the proton charge, Zy and Zz denote the charge num- 
ber of the heavy and light fragments, respectively, and den 
is the relative distance between the centers of charge of the 
two fragments at the scission point. This approximation for 
TKE has been frequently used for simplification, as demon- 
strated in Refs. [15, 20, 22, 57-59]. However, it neglects 
the dissipation and shell effects, among others, which can 
lead to an overestimated TKE compared to the experimen- 
tal data [15, 20, 22, 57-59]. The dissipation effect has been 
recently considered [60], thus the calculated TKE can bet- 
ter agree with the data. The TKE values of the ?40Pu fission 
fragments with different types of pairing forces were plotted 
as functions of the heavy fragment mass, as shown in Fig. 8. 
The open circles represent the calculated results for different 
7 values, whereas the solid circles indicate the experimental 
data from the thermal neutron-induced ?°9Pu fission experi- 
ments [61-63]. Considering the general trend, a qualitative 
dip is reproduced for all types of pairing forces at Ay = 120, 
as well as a peak at Ay = 134.Near the dip or peak, the TKE 
is larger for larger 7) values. For Ay > 144, the discrepancies 
are fairly small for different pairing forces. 


F. Fission yields 


Fig. 9 presents the mass and charge yields obtained with 
the code FELIX (version 2) [47] based on the TDGCM-GOA 
framework using the different types of pairing forces, which 
are compared with the experimental data. As a critical micro- 
scopic input of fission dynamic calculations, the mass tensor 
is calculated by the GCM or ATDHFB methods. In this calcu- 
lation, gv = 4 was used to determine the scission line. Gen- 
erally, the discrepancies between the calculated pre-neutron 
mass distribution and charge distributions obtained by us- 
ing the mass tensors by the GCM and ATDHFB methods are 
small. Furthermore, the mass and charge yields calculated by 
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Fig. 9. (color online) Calculated pre-neutron mass yields [panels (a) 
and (b)] and charge yields [panels (c) and (d)] are compared with 
the experimental data. Only the heavy fragments are displayed. In 
panels (a) and (c), the mass tensor from the GCM method is used, 
and in panels (b) and (d), the mass tensor is calculated using the AT- 
DHFB method. Data for the pre-neutron mass yields were obtained 
from Refs.[64, 65], and those of the charge yield were obtained from 
Ref. [66]. 


using the mixed-pairing force with 7 =0.5 combined with the 
mass tensor by the ATDHFB method demonstrated the best 
agreement with the experimental data. 

The impact of the different types of pairing forces on the 
mass and charge distributions is apparent. For the calculated 
results obtained by the ATDHFB mass tensor, the position of 
the peak was nearly constant for 7 = 0.0, 0.25, and 0.5, and 
moved toward the heavy fragment as 7 = 0.75 and 1.0. For 
the results with GCM mass tensor, the mass and charge distri- 
butions of the fission fragment shifted toward the more heavy 
fragment as 7) increased (panels (a) and (c)). Furthermore, the 
theoretical calculations obtained by the TDGCM with ATD- 
HFB mass tensors (panels (b) and (d)) demonstrated that the 
yields from the symmetric fission channel increased as n in- 
creased, which was related to the decrease in the height of the 
ridge as 77 increased, as shown in Fig. 2. 
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